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Abstract 

Signal-processing tools written in the Faust lan¬ 
guage are described. Developments in Faust libraries, 
oscillator.lib, filter.lib, and effect.lib since 
LAC-2008 are summarized. A good collection of sinu¬ 
soidal oscillators is included, as well as a large variety 
of digital filter structures, including means for specify¬ 
ing digital filters using analog coefficients (on the other 
side of a bilinear transform). Facilities for filter-bank 
design are described, including optional delay equal¬ 
ization for phase alignment in the filter-bank sum. 

Keywords 
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1 Introduction 

The Faust (Functional AUdio STream) language, 
developed at GRAME0p[2] , is well known for its 
compact specification of signal-processing block 
diagrams, and its compilation into efficient C-|—|- 
audio applications and plugins. Thanks to the 
use of architecture files that encapsulate platform- 
specific details, Faust applications can be conve¬ 
niently generated for a wide variety of host envi¬ 
ronments (Linux, Mac, Windows), and audio plu¬ 
gins can be generated for a wide variety of host ap¬ 
plications such as Pd and SuperCollider, to name 
just two PSIIS]. 

In the architecture subdirectory within the 
Faust distribution, there are presently seven 
. lib files containing various utility functions. 
Possibly the most commonly used of these is 
music.lib, which also imports math.lib. Three 
other . lib files pertain more specifically to signal 
processing utilities: 

*http://faust. grame.fr/ 


• oscillator.lib — signal sources 

• filter. lib — general-purpose digital filters 

• effect .lib — digital audio effects 

The remaining two Faust library files are 
maxmsp.lib—a Max/MSP compatibility library, 
and reduce.lib—enabling function application 
across a signal in time, such as maxn(n) = 
reduce (max,n) to compute the maximum ampli¬ 
tude of a signal. 

The directory examples/faust-stk/ addition¬ 
ally contains instrument. lib, providing com¬ 
mon utility functions for the Faust-STK collec¬ 
tion [6], such as envelope generators and table- 
lookup utilities. 

The libraries oscillator | filter | effect. lib 
were first discussed at LAC-08 [7]. This pa¬ 
per provides an overview of developments since 
then and up to Faust release version 0.9.46 
(Dec. 2011). 

2 Faust Library oscillator. lib 

The purpose of oscillator. lib is to provide 
reference implementations of various elementary 
waveform generators, such as sinusoidal, saw¬ 
tooth, and pulse-train, as well as other classic sig¬ 
nals such as pink-noise, etc. 

2.1 Sinusoid Generators 

All sinusoidal oscillators in oscillator. lib are 
invoked via the same API as osc(freq) (defined 
in music.lib), where freq is the desired osc¬ 
illation frequency in Hz. However, some provide 
two outputs instead of one when both “in-phase” 
and “quadrature” (sine and cosine) are available. 
All are filter-based. That is, they are implemented 
as lossless second-order filters driven by an im¬ 
pulse signal [1,0,0, ...], and they use no wave 



tablesll All algorithms have been previously pub¬ 
lished [SlElIiniEIli 

Presently, the following algorithms are imple¬ 
mented: 


oscb 

“biquad” two-pole filter section 
(impulse response) 

oscr 

2D vector rotation 
(second-order normalized ladder) 
provides sine and cosine outputs 

oscrs 

sine output of oscr 

oscrc 

cosine output of oscr 

oscs 

state variable osc., cosine output 
(modified coupled form resonator) 

oscw 

digital waveguide oscillator 

oscws 

sine output of oscw 

oscwc 

cosine output of oscw 


• oscs, based on the classic “state variable fil¬ 
ter,” m p. 530] and known as the “magic 
circle algorithm” in computer graphics, is 
quite fast, requiring only two multiplies and 
two additions per output sample. Its ampli¬ 
tude varies much less with frequency, and it 
too is good down to dc. There is no am¬ 
plitude drift over time, so this one can be 
used for very long signal durations. On the 
other hand, there is some dependence of osc¬ 
illation amplitude on frequency. At low fre¬ 
quencies, its two state variables are nearly in 
phase quadrature, but they become in-phase 
at SR/2. Thus, two outputs with approxi¬ 
mately 90-degrees relative phase at low fre¬ 
quencies could be brought out. The output 
that is brought out is the “cosine” choice. 


The relative merits of each oscillator type are 
summarized below. Note that all differences have 
to do with finite numerical precision effects and 
dynamic range variations under time-varying con¬ 
ditions. The best overall choice depends on the 
situation. 

• oscb, the impulsed direct-form biquacfl is the 
fastest computationally, requiring only one 
multiplication and two additions per sample 
of output. However, as is well known, the 
amplitude of oscillation varies strongly with 
frequency, and it becomes numerically poor 
toward freq=0 (“dc”). 

• oscr, the “2D vector rotation,” requires four 
multiplies and two additions per sample. Its 
amplitude is invariant with respect to fre¬ 
quency, and it is good all the way down to dc. 
Since its coefficients are numerically inexact 
roundings of s = sin(2*PI*freq/SR) and c 
= cos (2*PI*f req/SR), where SR denotes the 
sampling rate (defined in music.lib), there 
is long-term amplitude drift corresponding to 
the extent the identity s^-|-c^ = 1 is violated. 
This oscillator provides in-phase (cosine) and 
phase-quadrature (sine) outputs. 

^osc(freq) in music.lib uses a length 2^® wave table. 
The linearly interpolated variant osci(freq) adds linear 
interpolation. 

®https: //ccrma.Stanford.edu/~jos/pasp/- 
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• oscw, the second-order digital waveguide os¬ 
cillator, requires one multiply and three ad¬ 
ditions when frequency is constant, and an¬ 
other multiply when frequency is changing. 
Otherwise it has all of the good properties of 
oscr (except for internal dynamic range nor¬ 
malization), providing sine and cosine out¬ 
puts in exact phase quadrature, and no de¬ 
pendence of amplitude on frequency. How¬ 
ever, unlike oscr, oscw exhibits no amplitude 
drift while frequency is fixed. This is because 
it uses a “structurally lossless” algorithm de¬ 
rived by transformer coupling of normalized 
digital waveguides uni [8]Jj A negative point 
relative to oscr is that numerical difficulties 
may arise below 10 Hz or so, implying that 
oscw is not a good choice for LFOs. Inter¬ 
nally, the state variables of oscw require a 
larger dynamic range than those of oscr. It 
is likely that oscw would be the most eco¬ 
nomical choice for special-purpose VLSI. 

2.2 Virtual Analog Waveforms 

The following waveform generators are presently 
included, among others: 


imptrain(freq) 
squarewave(freq) 
sawtooth(freq) 
sawN(freq) 


periodic impulse train 
zero-mean square wave 
alias-suppressed sawtooth 
order N anti-aliased saw 


®https: //ccrma.Stanford.edu/~jos/pasp/- 
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The sawtooth and sawN algorithms are based 
on recently developed “Differentiated Polynomial 
Waveform” (DPW) methods for virtual analog 
waveform generation M M- The default case 
is sawtooth = saw2, where saw2 is a differenti¬ 
ated parabolic waveform (order 2). More gen¬ 
erally, sawN is based on a differentiated polyno¬ 
mial of order N. The higher the order, the less 
aliasing is incurred. Bandlimited square, triangle, 
and pulse-train are derived as linear filterings of 
bandlimited sawtooth in Faust releases beyond 
0.9.46. 

2.3 Noise Generation 

The basic white-noise generator, uniformly dis¬ 
tributed between —1 and 1, is noise, defined 
in music.lib. Based on that, oscillator.lib 
also defines pink_noise, also called “1// noise” 
implemented (approximately) as white noise 
through a three-pole, three-zero HR hlter that ap¬ 
proximates a 1// power response!! The third- 
order HR filter was designed using invfreqz in 
Octave (matlab). 

3 Faust Library filter, lib 

Filter-related utilities are provided in 
filter.lib. The principal functions defined 
appear in Fig.[Tl p. 01 and Fig.[2l p. [5l To save 
space, functions introduced at LAC-08 [7] are not 
repeated here (such as EKS string synthesizer 
elements, comb filters, cubic distortion overdrive, 
dc blockers, speaker bandpass. Crybaby wah 
pedal, etc.). The subsections below provide 
further discussion of various groups. The source 
is documented with comments and references so 
that anyone knowledgeable in basic digital filter 
theory [16] should be able to use it as a (terse) 
manual and starting point for further reading. 

3.1 Direct-Form Digital Filters 

The four direct-form digital filter structures have 
coefficients that appear in the filter transfer 
function]! The functions tfl(bO,bl,al) and 
tf2(b0,bl,b2,al,a2) specify first- and second- 
order (biquad) digital filter sections, respectively. 
Often larger filters are built by stringing first- 

®https: //ccrma.Stanford.edu/~jos/sasp/- 
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and second-order sections in series and/or par¬ 
allel. The Faust language makes this especially 
easy to do. 

The function iir(bcoeffs,acoeffs) allows 
specification of an arbitrary-order HR digital fil¬ 
ter in direct form. The arguments bcoeffs and 
acoeffs are each parallel signal banks that pro¬ 
vide the coefficients, and they may be thought 
of as “lists” of coefficients. The pattern match¬ 
ing facility in Faust allows recursive defini¬ 
tion in terms of such lists. As an example, 
tf2(b0,bl,b2,al,a2) can be alternatively spec¬ 
ified as iir((b0,bl,b2),(al,a2)). As usual in 
Faust, the specification is compact: 

iir(bv,av) = sub ~ fir(av) : fir(bv); 

where f ir(bv) specihes a general causal FIR dig¬ 
ital filter, with bv the list (“vector”) of FIR “tap” 
coefficients. It is given by 

fir(bv,x) 

= sum(i,count(bv),take(i+1,bv) * xOi); 
where count and take are defined in math.lib. 

3.2 Ladder and Lattice Digital Filters 

Ladder and lattice digital filters have superior nu¬ 
merical properties. Using the pattern-matching 
facility, it was possible to specify all four ma¬ 
jor types recursively in Faust. A particularly 
valuable case is the normalized ladder filter m 
iir_nl(bcoeffs,acoeffs), used as the basis for 
the super-robust biquad tf2snp(). While nor¬ 
malized biquads are straightforward to design 
{e.g., nlf 2() in filter. lib), the normalized lad¬ 
der filter realizes rational transfer functions of any 
order (any number of poles and zeros) in terms of 
a power-normalized ladder structure. For an in¬ 
troduction and pointers to references, see |8| and 
filter. lib|! 

3.3 Digital Filter Sections Specified as 
Analog Filter Sections 

It is convenient to be able to specify basic hlter 
section in terms of analog hlter coefficients, as op¬ 
posed to the usual digital-hlter coefficients. This 
is easy to do in Faust by including a built-in bilin¬ 
ear transform^ This makes use of the wonderful 

®https: //ccrma.Stanford.edu/~jos/pasp/- 
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Direct-Form Digital Filters 
fir(bcoeffs) 
iir(bcoeffs,acoeffs) 
tfl(bO,bl,al) 
tf2(b0,bl,b2,al,a2) 

(EH p.ED 

general FIR digital filter 
general IIR digital filter 

first-order direct-form digital filter = iir((bO,bl), (al)) 
iir((b0,bl,b2),(al,a2)) 

Lattice/Ladder Filters 
iir_lat2(bcoeffs,acoeffs) 
iir_kl(bcoeffs,acoeffs) 
iir.latl(bcoeffs,acoeffs) 
iir_nl(bcoeffs,acoeffs) 
tf2np(b0,bl,b2,al,a2) 
nlf2(f,r) 

(EH P-ED 

two-multiply lattice digital filter 

Kelly-Lochbaum ladder digital filter 
one-multiply lattice digital filter 
normalized ladder digital filter 

biquad based on stabilized 2nd-order normalized ladder 
second-order normalized ladder digital filter, special API 

Analog-Specified Filters 
tf2s(b2,bl,b0,al,a0,wl) 
tf2snp(b2,bl,b0,al,a0,wl) 
tf2sb(b2,bl,b0,al,aO,wl,wc) 
tflsb(bl,bO,aO,wl,wc) 

(EH p.ED 

tf2 specified via s-plane (analog) coefficients 
tf2s using a protected normalized ladder filter for tf2 
tf2s plus a mapping to bandpass in the digital domain 
same as tf2sb but for first-order filter sections 

IIR Low/High/Band-Pass 
lowpass(N,fc) 
highpass(N,fc) 
bandpass(Nh,f1,fu) 
bandstop(Nh,f1,fu) 
lowpass3e(fc) 
lowpass6e(fc) 
highpass6e(fc) 
bandpassl2e(f1,fu) 
bandpass6e(f1,fu) 

(E31 p. ED 

Ath-order Butterworth lowpass, —3 dB frequency at fc Hz 

Ath-order Butterworth highpass, —3 dB frequency at fc Hz 

Order 2*Nh Butterworth bandpass, —3 dB frequencies fl,fu Hz 

Order 2*Nh Butterworth bandstop filter, —3 dB gain at fl,fu Hz 

3rd-order elliptic lowpass, 60 dB stopband rejection, 0.2 dB passband rip. 
6th-order elliptic lowpass, 80 dB stopband rejection, 0.2 dB passband rip. 
highpass transformation of lowpass6e (w l/w) 

bandpass transformation of lowpass6e 
bandpass transformation of lowpass3e 

Shells, Peaking Eqnalizers 
low_shelf1(LO,fX,x) 
low.shelf1_1(GO,fX,x) 
low.shelf3(LO,fx,x) 
low_shelf5(LO,fx,x) 
low.shelf 
high.* 
peak_eq(Lfx,fx,B) 
peak_eq_cq(Lf X, f x, Q) 
peak_eq_rni (Lf X, f x, w) 

See Faust example parametric.eq.dsp 

Ist-order shelf, dc gain LO dB, crossover to unity gain at fx Hz 
dc gain GO (linear), crossover to unity gain at fx Hz 

3rd-order low shelf 

5th-order low shelf 

= low_shelf3; // default = third-order case 
same high-shelf cases as for low.shelf 

2nd-order “peaking equalizer”, peak level Lfx dB, width B Hz at fx Hz 
Constant-Q 2nd-order peaking equalizer section, Q = fx/B 

Regalia-Mitra 2nd-order peaking equalizer section, w ~ PI*B/SR 

Fractional Delay Lines 
fdelayN(maxdelay, delay) 
fdelayNa(maxdelay, delay) 

(EH p. ED 

Nth-order FIR Lagrange-interpolated delay line, N=l,2,3,4 

Nth-order HR allpass-interpolated delay line, N=l,2,3,4 


Figure 1: Functions defined in filter.lib since LAC-08. See the source code for full usage documenta¬ 
tion and literature references. 










Filter Banks 
nith_octave_analyzer (0, M, f top, N) 


mth_octave_analyzer6e (M, f top, N) 
mth_octave_f ilterbank(0,M,ftop,N) 
mth_octave_filterbank.alt 
mth_octave_spectral_level 

mth_octave_spectral_level6e 
spectral_level 
half_octave_analyzer(N) 
half _o ct ave_fiIt erbank(N) 
octave_filterbank(N) 
octave.analyzer(N) 
analyzer(0,Ifreqs) 
filterbank(D,lfreqs) 
filterbanki(0,Ifreqs) 


(Ezi p.m 

N-band octave filter-bank, M band-slices per octave, 

Butterworth band-split order 0 (not 0, must be an odd integer), 
N = total number of bands (including dc and Nyquist), 
ftop = highest band-split crossover frequency (e.g., 20 kHz) 
uses order 6 elliptic band-split filters 
mth_octave_analyzer followed by delay equalizer 
dc-inverted variant (cheaper for odd 0) 

spectrum analyzer using mth_octave_analyzer (5) , displays 
(in bar graphs) the average signal level in each spectral band 
order 6 elliptic crossovers 

= nith_octave_spectral_level(2,10000,20); // simplest 
= mth_octave_analyzer6e(2,10000,N) ; 

= mth_octave_f ilterbank5(2, lOOOOjN); 

= mth_octave_f ilterbankS(1,10000,N); 

= mth_octave_analyzer6e(l, 10000,N) ; 

general analyzer, order 0 Butterworth crossovers at listed freqs 
analyzer (0, If reqs) : delay equalizer (allpass-complementary) 
Inverted-dc variant 


Figure 2: Filter-bank functions defined in filter.lib. See the source code for full usage documentation 
and literature references. 


feature of Faust that if the coefficients are con¬ 
stant, all expressions will compile away to leave 
numerical digital-filter coefficients. On the other 
hand, if a slider-control, say, is providing an ana¬ 
log coefficient, the bilinear transform will be com¬ 
puted in real time (at the control rate) from the 
controller by the compiled result. Normally a one- 
pole smoother such as smooth(0.99) is used to 
interface the final computed coefficient into the 
filtering computation at the full sampling rate. 

In particular, tf 2s (b2a,bla,b0a, ala, aOa, wl) 
equals tf2(b0d,bld,b2d,aid,a2d) specified in 
the analog domain, where a last-letter ‘a’ means 
‘analog’, and ‘d’ means ‘digital’. (Note the oppo¬ 
site numbering of the coefficients, in conformance 
with typical notation.) Thus, the analog transfer 
function specified is 

N ^2,a + bi a s + 6o a 

H{s) = --. 

+ 0.1,a S + ao,a 

The parameter wl is the digital frequency ujd to 
which analog frequency Wa = 1 is mapped; it de¬ 
termines the frequency-scaling parameter of the 
bilinear transform. In lowpass or highpass filter 
design, the frequency mapping is applied to the 
cutoff frequency (—3 dB point). 


Butterworth filters are particularly easy to 
specify in analog form [THl [121 [S]) LJ because, 
for order N, all N zeros are at infinity and all 
N poles lie along a circle in the left-half s-plane. 
For example, the second-order Butterworth low- 
pass filter with its —3 dB frequency normalized 
to Wq = 1 is simply 


H{s) 


1 

s‘^ + V2s + l 


and can be specified as tf2s(0,0,1,sqrt (2), 1). 

3.3.1 Normalization and Stability 
Protection 

For extreme time-varying filtering applications, a 
practically useful variant named tf2snp is pro¬ 
vided that implements tf2s using a normalized 
ladder filter (for decoupling signal and coeffi¬ 
cient energy, ^3.2p together with stability pro¬ 
jection (easy to do in ladder/lattice digital fil¬ 
ters by simply clipping their reflection coefficients 
to the range (—1,1)). This is used in the most 
numerically robust Moog VCF implementation 
moog_vcf_2bn (effect.lib, 

https: //ccrma. Stanford. edu/~ jos/filter s/- 
Butterworth_Lowpass_Design.html 









The example vcf _wah_pedals. dsp in the 
Faust distribution provides a comparison of three 
Moog VCF implementations as well as the second- 
order Crybaby wah-pedal and a fourth-order wah- 
pedal based on the Moog VCF. 

3.3.2 Bandpass Mapped Biquad 

The function tf2sb(b2,bl,b0,al,a0,wl,wc) is 
a bandpass mapping of the basic analog-specihed 
biquad tf 2s. In addition to the frequency-scaling 
parameter wl (which gets set to half the desired 
passband width in radians per second), there is 
a desired center-frequency parameter wc (also in 
rad/s). Thus, tf2sb implements a fourth-order 
digital bandpass filter section specified by the co¬ 
efficients of a second-order analog lowpass proto¬ 
type section. Such sections can be combined in 
series for higher orders. The order of mappings is 

(1) frequency scaling (to set lowpass cutoff wl), 

(2) bandpass mapping to wc, then (3) the bilinear 
transform, with the usual scale parameter 2* SR, 
where SR denotes the sampling rate. The Faust 
implementation for this was based on algebra car¬ 
ried out in maxima. 

3.4 Butterworth Lo/Hi/Bandpass Filters 

Butterworth lowpass and highpass filters of any 
order can be dehned recursively in Faust thanks 
again to the pattern-matching facility in the lan¬ 
guage. The elliptic (Cauer) filterqlj are special- 
cased because the pole locations are computed 
using the elliptic rational function, which is not 
available in typical computer math libraries. Such 
a function could of course be supplied as a foreign 
function in Faust. 

3.5 Shelf and Equalizer Sections 

The low/high shelf and peaking equalizer sections 
implemented in filter. lib are described further 
in filter.lib and in [T6]0 

3.6 Lagrange/Thiran-Interpolated 
Fractional Delay Lines 

Delay lines interpolated using higher-order FIR 
Lagrange interpolation are all used as follows: 

fdelayN(maxdelay, delay, inputsignal) 

^'http://en.Wikipedia.org/wiki/Elliptic_fliter 
^https: //ccrma.Stanford.edu/~jos/fliters/- 
Low_High_Shelf-Filters .html 


where N=l,2,3, or 4 is the order of the Lagrange 
interpolation polynomial. Note that this API fol¬ 
lows that of f delay in music. lib. The requested 
delay should not be less than {N — l)/2 because 
the interpolating polynomial needs to be able to 
“reach” that far into the “past” when interpolat¬ 
ing. 

Delay lines interpolated using higher-order HR 
allpass Thiran interpolation are all invoked as 

fdelayNafmaxdelay, delay, inputsignal) 

where N=l,2,3, or 4 is the order of the allpass in¬ 
terpolation filter. In this case, it is recommended 
that the requested delay be at least N — 1/2 be¬ 
cause an Vth-order allpass provides a delay of N 
samples as its coefficients approach zero. Note 
that delay arguments that are too small can pro¬ 
duce an unstable allpass filter. For rapid delay 
modulations, Lagrange (FIR) interpolation is gen¬ 
erally preferred. However, allpass interpolation 
introduces no gain distortion and may therefore 
be preferred in nearly lossless feedback loops. 

3.7 Filter Banks 

A filter bank splits its input signal into a bank 
of parallel signals, one for each spectral band. 
If the bandpass hlters used to create the chan¬ 
nel signals are carefully designed, one may sum 
the channel signals to get back the original input 
signal (possibly scaled and/or delayed). In this 
case, the filter bank is said to be a Perfect Re¬ 
construction (PR) filter bank |20j . However, for 
purposes of spectrum analysis, in which only the 
channel signal powers are displayed, the PR con¬ 
dition is overkill. Therefore, the filter banks im¬ 
plemented in filter.lib are divided into “an¬ 
alyzers”, which do not have the PR property, 
and “filter banks” which are “allpass complemen¬ 
tary”. Allpass-complementary filter banks are 
reasonable choices for “graphic equalizer” appli¬ 
cations. An allpass-complementary filter bank is 
PR when the allpass reduces to a pure delay and 
possible scaling. In this terminology, the hlter 
banks in filter.lib are implemented as analyz¬ 
ers in cascade with delay equalizers that convert 
the analyzer to an allpass-complementary hlter 
bank. Spectrum analyzer outputs should at least 
be nearly “power complementary”, i.e., the power 
spectra of the individual bands should at least ap¬ 
proximately sum to the original power spectrum. 



The typical filter bank or analyzer is con¬ 
structed as a dyadic filter bank, meaning that 
it consists of a sequence of band-splits, forming 
a binary tree of lowpass/highpass filter sections. 
Since audio applications are presumed, only the 
lower band is split when going from one stage to 
the next. 

In the Faust distribution, both filter banks and 
spectrum analyzers are illustrated in the example 
graphic_eq.dsp. See also spectral_level.dsp 
which is a standalone spectrum analyzer (nice as 
a standalone JACK app). 

The example gate_compressor. dsp included 
with the Faust distribution exercises the gate 
and compression utilities. 

Space limitations preclude further discussion 
here. Please see comments in filter. lib for fur¬ 
ther usage details. 

4 Faust Library effect, lib 

The modules in effect.lib classify as “digital 
audio effects”. In general, they tend to be special- 
purpose filters, frequently nonlinear and/or time 
varying. 

4.1 Moog Voltage Controlled Filters 

New since the analog-form Moog VCF [7] is the 
implementation moog_vcf_2b of the ideal Moog 
VCF transfer function factored into second-order 
sections. As a result, its static frequency re¬ 
sponse is more accurate than moogwcf which 
suffers from an unwanted one-sample delay in 
its feedback path. On the downside, its coef¬ 
ficient formulas are more complex when one or 
both parameters are varied. The res parameter 
of moog.vcf _2b [n] is the fourth root of that in 
moog_vcf , so, as the sampling rate approaches in¬ 
finity, moog_vcf (res,fr) becomes equivalent to 
moog_vcf_2b[n](res“4,fr) (when res and fr 
are constant). 

4.2 Artificial Reverberation 

The reverberation modules in effect. lib are de¬ 
scribed in |8]. Of special note is the high-quality 
reverberator called zita_revl, ported to Faust 
from the C-|-+ source of zita-revl written by Fons 
Adriaensen0 It combines Schroeder allpass and 
FDN reverberation techniques [8]P^ 

^^http://kokkinizita.linuxaudio.org/linuxaudio/- 
zita-revl-doc/quickguide.html 

^"^https: //ccrma. stanford.edu/~jos/pasp/Zita_Revl .html 


5 Conclusion 

Developments since LAC-08 for Faust libraries 
oscillator I filter I effect. lib were outlined. 
The overall goal is to accumulate reference imple¬ 
mentations of commonly used algorithms in mu¬ 
sic /audio signal processing, with a general prefer¬ 
ence for expressive parametric algorithms yielding 
the highest performed sound quality per unit of 
computation. 
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